Levy Flights from a Continuous-Time Process 
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The Levy-flight dynamics can stem from simple random 
walks in a system whose operational time (number of steps n) 
typically grows superlinearly with physical time t. Thus, this 
processes is a kind of continuous-time random walks (CTRW), 
dual to usual Scher-Montroll model, in which n grows sublin- 
early with t. The models in which Levy-flights emerge due to 
a temporal subordination let easily discuss the response of a 
random walker to a weak outer force, which is shown to be 
nonlinear. On the other hand, the relaxation of en ensemble 
of such walkers in a harmonic potential follows a simple expo- 
nential pattern and leads to a normal Boltzmann distribution. 
The mixed models, describing normal CTRW in superlinear 
operational time and Levy-flights under the operational time 
of subdiffusive CTRW lead to paradoxical diffusive behavior, 
similar to the one found in transport on polymer chains. The 
relaxation to the Boltzmann distribution in such models is 
slow and asymptotically follows a power-law. 

PACS No.: 02.50.-r, 05.40.Fb, 82.20.M 



I. INTRODUCTION 

Random walk processes leading to subdiffusive or su- 
perdiffusive behavior are adequate for describing vari- 
ous physical situations. Thus, the continuous-time ran- 
dom walk (CTRW) model of Scher and Montroll § 
was a milestone in understanding of photoconductivity 
in strongly disordered and glassy semiconductors, while 
the Levy-flight models Q are adequate for description 
of transport in heterogeneous catalysis (|] , self-diffusion 
in micelle systems 0], reactions and transport in poly- 
mer systems under conformational motion |fj], transport 
processes in heterogeneous rocks ||, and for description 
of behavior of dynamical systems 0. The closely re- 
lated models appear in description of economic time se- 
ries H . The Levy-related statistics were observed in hy- 
drodynamic transport Jl , and in the motion of gold nan- 
oclusters on graphite [H]]. The mixed models were pro- 
posed , in which the slow temporal evolution (described 
by Scher-Montroll CTRW) is combined with the possi- 
bility of Levy-jumps, so that in general both sub- or su- 
perdiffusive behavior can arise lit 

The continuous-time random walks (CTRW) first in- 
troduced by Montroll and Weiss [jl2| correspond to a 
stochastic model in which steps of a simple random walk 
take place at times ti, following some random process 
with non- negative increments: Tj = ti — ti—i > 0. In 
a mathematical language one says that CTRW is a pro- 
cess subordinated to random walks under the operational 



time defined by the process {ti}. It is typically thought 
that a CTRW-scheme alone can not describe any su- 
perdiffusive process, so that the introduction of very long 
jumps is an inevitable part of building a model leading 
to superdiffusive behavior. 

Let us first discuss a typical CTRW approach. Let us 
consider a one-dimensional situation under which a par- 
ticle from time to time makes a jump to a neighboring 
lattice site separated from the initial one by a distance 
a. The time r between the two jumps is distributed ac- 
cording to some waiting-time distribution, represented 
by the probability density function (PDF) p(r). If the 
mean waiting time f exists, the particle's behavior is 
diffusive, with diffusion coefficient D = a 2 /2f. If the 
corresponding moment diverges, the particle's behavior 
becomes subdiffusive, with (r 2 (t)) oc t a , with a < 1 de- 
pending on the PDF p(r). The subdiffusive behavior is 
indicated by vanishing of the diffusion coefficient D. It 
seems impossible to obtain within this scheme any type of 
a superdiffusive behavior unless one allows for infinitely 
long jumps with (a 2 ) — * oo. The superdiffusive behavior 
is indicated by divergence of the diffusion coefficient D. 
If (a 2 ) stays finite this can be the case only if f van- 
ishes. Since r > and f = J °° rp(r)dr, vanishing of the 
mean waiting time means that p(r) = 5{r), a marginal, 
degenerate situation. 

On the other hand the consideration presented above 
shows only that the waiting-time distribution is not an 
adequate tool for description of superdiffusive CTRW. 
In what follows we show that superdiffusive CTRW with 
bounded step lengths are just as possible as the subdiffu- 
sive ones. Our considerations will be rather formal and 
do not follow from any particular physical model. On the 
other hand, the fact that Levy-flights can stem from a 
process subordinated to simple random walks has many 
important implications. Thus, as we proceed to show, 
the fast dynamics of a free process can coexist in such 
models with simple exponential relaxation to a normal 
Boltzmann equilibrium distribution, if the behavior of an 
ensemble of random walkers under restoring force is con- 
sidered. This shows that the relation between Levy dy- 
namics and the nonextensive thermodynamics described 
by nonclassical entropy functions is much looser than typ- 
ically assumed. 

The combinations of the superdiffusive Levy-flights 
with the typical CTRW operational time leads to para- 
doxical diffusion behavior, having some parallels in trans- 
port on polymer chains. Moreover, the existence of a sub- 
ordination model leading to Levy flights can be useful in 
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understanding of statistical implications of the processes 
described by fractional generalizations of diffusion and 
Fokker-Planck equations [p~3| — . 

The article is organized as follows: In Sec. 2 we discuss 
general properties of subordinated random processes. In 
Sec. 3 and 4 the processes subordinated to symmetric and 
asymmetric random walks are considered, these leading 
to symmetric and asymmetric Levy-flights. The dualism 
between the Levy-flights and the Scher-Montroll CTRW 
is discussed in Sec. 5. Sections 6 and 7 discuse the models 
leading to paradoxical diffusion behavior. The relaxation 
to equilibrium is considered in Sec. 8. 

II. THE SUBORDINATION OF RANDOM 
PROCESSES 

As already mentioned, a Scher-Montroll CTRW pro- 
cess is a simple random walk whose steps take place at 
times tj governed by a random process with nonnegative 
independent increments, so that 

P(x > t) = ^P RW (x,n) Pn (t), (1) 

n 

where Puw(x,n) is a probability distribution to find a 
random walker at point x after n steps (i.e. the binomial 
distribution), andp„(t) is the probability to make exactly 
n steps up to time t. For both t and n large, when the 
binomial distribution can be approximated by a Gaussian 
one, and when the corresponding sum can be changed to 
an integral, Eq.([l]) reads: 

f°° 1 / x 2 \ 

P(x,t)~ / -7=exp -— )p t (n,t)dn. (2) 

In a classical Scher-Montroll CTRW pt(n, t) corresponds 
to a random process in which n typically grows sublin- 
early in t. Thus, the overall process is subdiffusive. 

Note that a description of CTRW-process given by 
Eq.( B) i s an example of subordination, see Sec. X.7 of 
Ref . |lq| : If {X(T)} is a Markov process with continuous 
transition probabilities and {T(t)} a process with non- 
negative independent increments, then {X(T(t))} is said 
to subordinate to {X(t)} using the operational time T. 
In this case 

poo 

P(x,t) = / P x (x,T) PT (T,t)dT. (3) 
Jo 

In what follows we call the integral transform, Eq.(||) a 
subordination transformation, changing from time scale 
i to a time-scale T. For example, in the Scher-Montroll 
case the operational time of a system is given by the 
number of steps of the RW, and is a random function of 
the physical time t whose typical value grows sublinearly 
in t. 

The operational time can also grow superlinearly with 
t. Such a process can not be described by a waiting-time 



distribution, and needs a complimentary description. Let 
us consider a random process, where the density of events 
fluctuates strongly. Let us subdivide the time axis into 
intervals of duration At and let us consider the num- 
ber n of jumping events within each interval. The value 
p = n/ At defines the density of jump events. Now, if the 
mean density of events exists, its inverse gives us exactly 
the mean waiting time of a jump, and a process described 
by a finite density of events is a normal diffusive one. The 
divergence of a mean waiting time (like in Scher-Montroll 
CTRW) correspond to vanishing density. On the other 
hand, if one considers a strongly fluctuating density p(t) 
whose first moment diverges, the mean waiting time van- 
ishes and a process that subordinates a random walk pro- 
cess under such operational time can be superdiffusive. 
At longer times, the distribution of the number of events 
tends to one of the Levy-stable laws: the typical number 
of events can grow superlinearly in time. A simple ex- 
ample of such process was already known to Feller, see 
Chap. X.7 of Ref. He considers a process subor- 

dinated to simple random walks under the operational 
time governed by a fully asymmetric Levy stable law of 
index 1/2. The corresponding PDF at time t is given by 

p(M) = i 72b exp (-^)7^ exp (^) dr ' 



ir(t 2 + x 2 )' V 

i.e. is a Cauchy Levy-flight. 

Let us now discuss a simple analogy describing the 
relation between the Scher-Montroll CTRW and Levy- 
flights. This analogy makes clear many of the findings 
we are going to discuss below. Imagine a physical clock 
producing ticks following with frequency 1, which gov- 
ern the behavior of a random walker. Imagine a switch 
situated at 0, so that returning to the origin, the walker 
can trigger some physical process (the analogy with the 
Glarum model of relaxation, Ref. |17|| , is evident!). The 
times between the subsequent returns are distributed ac- 
cording to a fully asymmetric Levy stable law of index 
1/2 used in a previous example. Imagine now another 
random walker performing its motion (a step per physi- 
cal unit time) independently from the first one. Imagine 
a movie camera, taking frame-per-frame pictures of the 
positions of this second random walker at the moments 
when the first walker is at the origin and thus triggers the 
switch. Watching the movie taken by the camera, we im- 
mediately recognize that the second walker performs the 
Cauchy Levy-flights. Imagine, that a clock is posed in a 
frame and also filmed. In this case its image will show 
exactly the operational time of the system; the specta- 
tor's watch measures the physical time. Imagine an oppo- 
site situation: the first walker triggers the motion of the 
second one, and the camera is triggered by the physical 
clock, as a normal movie camera is. The process we rec- 
ognize at the film is then the Scher-Montroll CTRW. We 
can take a Scher-Montroll movie also using another trick 
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(which cannot be performed in a real time, but needs 
a record of return times). Let us take a record of sub- 
sequent return times of a first random walker (numbers 
m, ri2, ■■■) and trigger our camera in such a way that 
it makes n\ frames during the first second, frames 
during the second one, etc. If we film a normal random 
walker with a camera prepared in such a way, the movie 
will show us the Montroll- Weiss CTRW. An image of the 
physical clock will again show the operational time of the 
system, and again, looking at his watch, the spectator 
can measure the physical time between two events. 

Let us use our camera triggered by returns of a ran- 
dom walker to film other processes taking place in the 
outer world. The film which is watched afterwards under 
constant speed shows us a possible world: The causality 
relations and thermodynamical time arrow are those of 
our usual world. On the other hand, a movie of a world 
undergoing continuous evolution, in which " natura non 
facit saltus" holds, will show us a revolutionary world 
of "great leaps" and abrupt changes (but following the 
same logics of development). The second camera (fed by 
a prescribed n-sequence) will show us the world of almost 
full stagnation seldomly interrupted by a bounded, local 
movement, a world developing in a slow time of old Asi- 
atic despoty. We shall keep this analogy in mind when 
discussing the physical implications of subordination. 

Let us consider a system which evolves according to a 
Markovian dynamics and whose state tends to a normal 
Boltzmann equilibrium under relaxation. In a system 
under action of outer forces, the transition probabilities 
between the states of the system (sites i between which 
the random walk takes place) which are characterized by 
their energies J3j, are not independent. They are con- 
nected through the corresponding Boltzmann-factors, so 
that in equilibrium during any period of time At the 
mean numbers of forwards and of backwards jumps be- 
tween any two sites i and j fulfill the condition 



nij(At)/ nji (At) = exp - E 3 ) /kT] . 



(5) 



where k is the Boltzmann constant and T is the sys- 
tem's temperature. The condition Eq.(||) guarantees de- 
tailed balance in equilibrium, independentl of what the 
real dynamics of a system is. For simple RWs, where 
only transitions between the neighboring states are al- 
lowed, the corresponding transition rates with respect to 
the operational time of the system can be introduced. 
For a random a walker moving under the influence of a 
weak constant force F the probabilities of the forward 
and backward jumps per unit time w+ and W- are con- 
nected through w + /w- = exp(Fa/kT). The Markovian 
nature of RW then leads to the fact that the values of 
u>+ and W- do not depend on whether the system is 
in equilibrium or not. For F small one can take, say, 
w + — wq(1 + Fa/kT) and w_ = Wq(1 — Fa/kT) with 
w = 1/2t. 

Note that subordination, describing a transition from 
a physical time to an operational time of the system, 



does not change its equilibrium properties. Such sub- 
ordination can be considered as random modulation of 
the transition rate wq by some independent process (say 
closing and opening the channels), and is fully irrelevant 
for thermodynamics (i.e. thermostatics) of the system. 
On the other hand, it strongly influences its kinetics, so 
that a question can be posed, what kinds of kinetics are 
compatible with the relaxation to a normal Boltzmann 
distribution under arbitrary subordination transforma- 
tion of time. We address this question in Sec. 8, after the 
free diffusion properties of superdiffusive CTRW will be 
discussed. 



III. SYMMETRIC LEVY FLIGHTS FROM CTRW 

Let us first concentrate on the symmetric random walk 
case. Let us consider a random process in which the num- 
ber of events per given time is unbounded and follows, 
for example, a power-law distribution, p n (t) oc tn~ l ~ a 
with < a < 1 (this corresponds to the typical number 
of events scaling asn« t 1 ^"). Let us find the asymptotic 
behavior of P{x, i) for t large. Since the jumps during dif- 
ferent intervals are uncorrelated, the PDF of n for longer 
times converges to a fully asymmetric Levy-stable law 



p(n,t) ~ t _1 / a L(n/t 1/a ;a,7) 



(6) 



with the asymmetry parameter 7 = — a (here the val- 
ues of 7 = ±a correspond to the strongly asymmetric 
PDF that vanish identically for large positive (negative) 
x values , while 7 = corresponds to symmetric distri- 
butions; the notation in one of Ref. fll6fl ). Note that the 
Fourier-transforms of Levy-stable laws are known: up to 
the translation P(k, t) is equal to 



/(«) = exp - e 



1-7/2 



(7) 



(for 0<a<2, a^l). The PDF is a real function, thus 
f( K ) = /*( — The corresponding function is analytical 
everywhere except for k — 0, so that the PDF is given 
by " 



L(x; a, 7) = —Re 



(8) 



From Eq.^| the series expansions for L(y;a,j) follow, see 
Sec. XVII. 6 of Ref. Jig . In the case a < 1 one can 
move the path of integration to the negative imaginary 
axis (since the integrand tends to zero when Im£ — > — 00 
due to the dominance of the linear term), which allows 
then for elementary integration after Taylor-expansion 
of exp(A£ a ). For 1 < a < 2 this dominance is no more 
the case, but the integrand still vanishes for lm( — > — 00 
in the case of symmetric distributions, while (— i |C|) Q = 



-00 for 1 < a < 2. Thus the 



|C| a (cos %a 

series which represents Levy distributions for < a < 1 
and symmetric Levy distribution also for 1 < a < 2 
reads: 
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—ah 



(9) 

In general the Levy-stable laws for 1 < a < 2 are given 
by another expansion, 

ti \ 1 1 \fc r ( 1 + k / a ) ■ ( kn i \\ -k 

L(y;a, 1 ) = —J2(-V k fc , ' ' an ( — (7- a) J g 

y fc=i ' v 7 

(10) 

which also holds for asymmetric laws. 

One can easily obtain the form of x-distributions by 
immediate integration: using Eq.(|^) and a scaling form 
of a Levy-distribution 



p(x,t) 



] / x 2 \ . n . dn 

exp — — L( — a, — a 1^-7—. 



'27m 



(11) 



Using Eq. (^) and performing a term- by- term integration, 
we arrive to the series of integrals of the form 



WO 



1 (( 



2tt V 2 



1 c 2 



2tt£ 



r(/i-i/2). 



(12) 



For integral of the fc-th term in Eq.(||) we have /j, = 
1 + ak. Let us concentrate first on the case < 
a < 1. Using well-known relations for F-function: 
L(z + 1) = zT(z) (Eq. (6.1.15) of Ref. g|) and L(2z) = 
(27r)- 1 / 2 2 2z - 1 / 2 L(z)r(z + 1/2) (Eq. (6.1.18) of g§) we 
get 



fe=i 



(13) 



which represents a series expansion for a symmetric 
Levy-stable law of index 2a, Eq.(^J), for the scaled 
variable C/V2- This corresponds to a form p(x,t) = 
t- x l 2a L(xl^2t 2a ; 2a, 0) of the x-distribution. 

We note that taking Fourier-transform of the both 
parts of for symmetric RWs, 



L(ax, 2a, 0) = 



- 2 \ 



, exp ( I L(n; a, —a)dn, 

'2im V 2 " / 

(14) 



where a is an unimportant scaling factor, we get: 

/>oo 

exp(— A I A; I a t)= / exp (— fc 2 n) L(n; a, —a)dn, (15) 
Jo 



which holds for any real k (i.e. for any positive k 2 ), 
where A is a number factor. This gives us a general 
expression for a Laplace-transform of an asymmetric 
Levy-distribution with a < 1: A Laplace-transform of 
L(n,a,—a) is exp(— A\k\ t). From this fact an impor- 
tant result follows: 

/>oo 

L(ax; a/3, 0) = / n - 1/f3 L(x/n 1/l3 ; f3, 0)L{n; a, -a)dn : 
Jo 

(16) 

A Levy distribution with index a[3 is subordinated to 
a Levy-distribution with index j3 < a under the opera- 
tional time given by an asymmetric Levy law of index 
a < 1. To see this, consider the characteristic functions 
of both sides of Eq.© and use Eq.(||). 



exp(-^4 |fe| 



e-^ n L{n;a,-a)dn, (17) 



see Sec.X.7 of Ref. pi. Eq. (|15|) corresponds to a 
special case of [3 = 2 of Eq.Ql7j). The distribu- 
tions L(n; a, —a) thus coincide with inverse Laplace 
transforms of stretched-exponentials. For example for 
L{n; 1/2, —1/2) one readily gets: 



p(n,t)=£- 1 {expt-iu 1 / 2 )} 



t 



2^n i / 2 



ex P ~1- 



An) ' 
(18) 



which differs only by a scale for the time-unit from a 
distribution used in the example Eq.(0). 



IV. ASYMMETRIC LEVY-FLIGHTS 

Imagine a random walker moving under the influence 
of a weak constant force F. Such force introduces an 
asymmetry into the walker's motion, since the proba- 
bilities of the forward and backward jumps, w+ and 
W— are now weighed with the corresponding Boltzmann- 
factors, iv+/w- = exp(FafkT). For F small one can 
take w + = 1/2 + Fa/2rkT and w + = 1/2 + Fa/2rkT. 
For t large such random walks lead to the Gaussian dis- 
tribution of the particles' positions 



PRw(x,t) 



1 



'27m 



exp 



(x 



2n 



(19) 



whose center moves with a constant velocity v = fiF = 
Fa 2 /2rkT. Note that our RW fulfil the Einstein's rela- 
tion between the mobility \x and diffusion coefficient D: 
fi = D/kT. The PDF of a random process which subor- 
dinates biased RW under an operational time following 
the asymmetric Levy-law is given by: 



P(x,t) 



1 



'27m 



exp 



[x — vnY 
2n 



n dn 
xL{ . , ; a, —a . , . 



(20) 
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Using the series expansion, Eq.(|9j) and performing the 
term-by-term integration leads to the series of the inte- 
grals of the type: 



2exp(Cw) / C 



2tt 



2 x l/4-/i/2 



(Cw) (21) 



for to ^ 0. For integral of the A:-th term in Eq.(||) we 
again have /i = 1 + afc. Let us concentrate first on the 
case < a < 1. For small, u cancels (see the expan- 
sion 9.6.9 of Ref. @, K v {z) ~ ±r(z/)(±,z)- I/ (i/ > 0), 
note that -ftT_„(z) = K u {z)), so that the corresponding 
distribution tends to be a function of C, only, it coin- 
cides with one for ui = 0, Eq.jHj) so that a symmetric 
Levy-stable law of index 2a, Eq.Q emerges. On the 
other hand, for » / and x large the overall distribu- 
tions follow from the expansion of K for large values of 
the argument which reads: K v (z) ~ y/j; e ~ z (Eq. 9.7.2 
of Ref. 1 18 1). The corresponding integral then tends to 
■i (C/oj) so that the corresponding PDF reproduces 
the PDF of the density of events (up to reseating) . This 
last form is also the asymptotic from corresponding to 
the behavior of Eq.(p0|) for large t. 

Hence, the distribution P(x, t) tends to a fully asym- 
metric one of index a for x and t large. In this case 
the distribution shows scaling with a scaling parameter 
£ = x/(vt) a . We see that in this case the motion under 
the influence of a constant force is superdiffusive, so that 
x ~ (Ft) 1 / 01 , and its dependence on the outer force is 
nonlinear. Thus, the model shows a behavior that differs 
considerably from a linear-response assumption of Refs. 
p] , ^9|j20| . This absence of a linear response regime is 
parallel to the CTRW-findings @ (see Ref. fl| for a re- 
view) and mirrors the fact that only for normal diffusion 
a sweep with constant velocity and a drift under a con- 
stant force result in the same pattern of motion, see Ref. 
& 

The case a = 1/2 again results in a closed expression: 
I ( (x — vri) 



P(x,t) 




exp 



v27m J 
vt 



exp 



e_ 

2n 



2n 
dn 



(22) 



1 



7T y/x 2 + t 2 



e vx K Y [yjv 2 (x 2 + t 2 ) 



(2.3.16.1 of Ref. ]2~5[|). For v,x and t small, the corre- 
sponding distribution tends to a Cauchy-law. On the 
other hand, for t large we can take approximately: 



P(x,t) 



Vvt 



2tt ( x 2 + 1 2 ) 3/4 



(23) 



The second moment of this distribution diverges, but 
the position of the maximum of P(x,t), determining 



the typical particle position at time t, tends to grow as 
x max = |i 2 for t large. Thus, the typical behavior of x(t) 
under constant force is superlinear. 

Note that in the case 1 < a < 2 the distribution of 
the particle's displacement for the case v = will tend 
to a Gaussian, but in the case v > it still tends to a 
fully asymmetric Levy one. On the other hand, in this 
case the distribution of the particle's position possesses 
the first moment which grows linearly with time, thus 
the situation under a > 1 shows the linear response be- 
havior. Since the second moment of the distribution is 
absent, the fluctuations are strong, and the width of such 
distribution is of the order of the typical value of x itself. 



V. THE DUALISM BETWEEN THE SUB- AND 
THE SUPERDIFFUSIVE CTRW 

There exists a clear dualism between the normal, subd- 
iffusive CTRW and a superdiffusive one. The correspond- 
ing concepts are illustrated in discrete time by Fig.l, 
where we return to a situation discussed in Sec. 2. Imag- 
ine a clock producing ticks following with frequency 1, 
marking the physical time of a system. Imagine a system 
which is triggered not by each tick of a physical clock, but 
follows some waiting-time distribution, ip{n). This means 
that after our random walker has jumped, the next jump 
will take place after n ticks of a clock, where the num- 
ber n is chosen according to a power-law distribution, say 
ip{n) oc (n+l)~ 1+7 . The number n fluctuates strongly, so 
that the sequence of jumps (corresponding to a randomly 
decimated sequence of ticks) shows lacunae of different 
duration. Fig. la) shows a realization of such a sequence 
for the case 7 = 0.75. The lacuna starting in the middle 
of Fig. la) at t — 54 ends at t = 161. The mean num- 
ber of jumps during the time t grows sublinearly with t, 
namely as i 3 / 4 . Let us denote the corresponding subordi- 
nation transformation as time-expanding transformation 
(TET) of index 7. According to the procedure described 
above, the corresponding sequence does not have any in- 
tervals where the density of events is larger than one. 
The process subordinated to random walks under such 
operational time (normal CTRW) is subdiffusive. 

Let us now consider the sequence of jumps of a walker 
as ticks marking relevant time epochs of a system (i.e. 
associate each jump with a tick of a physical clock). 
From this point of view, the ticks of initial clocks fol- 
low extremely inhomogeneously, so that the number of 
such ticks within a physical time unit varies according 
to p(n) oc (n + 1)~ 1+7 . Fig. lb) illustrates this situation: 
Here we took 100 jumps from the realization shown in 
Fig. la) and rescaled each of the corresponding time inter- 
vals to the unit length. The ticks of initial clock (shown 
as bars) follow inhomogeneously and show the intervals 
of high concentration (but no lacunae). The number 
of such events grows superlinearly in time. The corre- 
sponding subordination transformation will be called a 
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" time-squeezing transformation" (TST) of index 7. The 
process subordinated to random walks under such oper- 
ational time is superdiffusive and corresponds to Levy- 
flights. Note that both TST and TET are the probability 
distributions P{n, t) of the operational time n for a given 
physical time t, i.e. are positive, integrable functions of 
n. 



a) 










operational 
















physical 




physical 



operational 



1 1 1 1 1 1 

20 40 60 80 100 



Fig. 1. This figure illustrates the notion of the operational 
time, a) The operational-time leading to the Schcr-Montroll 
CTRW. The "ticks" triggering the motion of the random 
walker (shown as bars above the horizontal line) are taken 
from the set of the "ticks" of a physical clock according to 
a waiting-time distribution ip{ri) oc (n + 1) — 1+7 (shown as 
bars below the line). The process shows long lacunae, see text 
for details, b) We now declare each interval between the two 
moves of the random walker for a new time unit. The initial 
ticks of the physical clock follow now extremely inhomoge- 
neously and show the intervals of high condensation. This 
kind of operational time leads to Levy-flights. 



Let us now combine the two processes. For example, 
let us first generate the superlinear sequence using algo- 
rithm described above (with p(n) oc (n + 1)~ 1+7 ) and 
then decimate it randomly according to the waiting-time 
distribution ip(n) oc (n + 1)~ 1+7 . The typical number 
of events during time interval t grows in this case lin- 
early with t, but the corresponding sequence of events is 
extremely inhomogeneous, showing both lacunae and ac- 
cumulation intervals on all scales. This process is shown 
in a bar-code-like picture in Fig. 2. It will be discussed 
in more detail in Sections 6 and 7. We can also pro- 
ceed other way around, and apply the transformations in 
the opposite way, namely, first generating a sublinearly 
growing, lacunary operational time and then filling the 
lacunae according to a Levy distribution. As we proceed 
to show, these two ways of constructing the event-time 
sets are not equivalent. The process subordinated to RW 
under such inhomogeneous operational time is a kind of 
a continuous-time Levy- flight, and not a normal RW. 

The example discussed above shows that transforma- 
tions leading to sub- or superlinear operational time be- 
havior (dual to each other in the sense described above) 
are not inverse of each other. Let us discuss a possibility 
of a subordination transformation transforming a Levy- 
stable distribution of index (3 (for example, a Gaussian 
distribution) into one with ones with index 7, in the sense 
that 



L(ax;j,0) = n - 1/,3 L(x/n 1/l3 ; (3, 0)S(n, t)dn, (24) 
Jo 

where S(n, t) is supposed to be a probability distribution 
of the number of steps n done up to time t. Taking 
Fourier-transform of both parts of Eq. (|24|) and changing 
to a variable u = Ik] 13 we get: 

/>oo 

exp(-A\u\ a ) = / e- un S{n,t)dn (25) 
J 

with a — 7//?. From Eq.(|25|) it follows that S(n,t) are 
the inverse Laplace transforms of stretched-exponentials 
exp(— Au a ). Note that according to the Bernstein's the- 
orem, a function f{x) is a Laplace-transform of a proba- 
bility distribution if and only if it is completely monotone 
(i.e. it is infinitely differentiable and (— l) n f( n \x) > 
for all derivatives / ( -™- ) ) and /(0) = 1. The last condi- 
tion is always fulfilled. Note that according to Criterion 
2 discussed on p. 441 of vol.IIII of Ref. p6[ a function 
f(x) = e~^' l - x ' > is a completely monotone function if and 
only if ip is a positive function with a completely mono- 
tone derivative. In our case ip( x ) = Au a . For < a < 1 
one has: g(x) = ip'(x) = Auu a ~ l > 0, and the higher 
derivatives (defined on the interval < x < 00) are: 
g'{x) = Aa(a - l)u a - 2 < 0, g"(x) = Aa{a - I) (a - 
2)u a - 3 > 0, g"'{x) = Aa(a - I) (a - 2)(a - 3)u a ~ 3 < 0, 
etc., so that (— l) n g^ n \x) > 0, and thus the function g is 
completely monotone. Thus S(n, t) is a probability dis- 
tribution (namely the one which we have found above by 
explicit calculation). On the other hand, for a > 1 the 
function g{x) is not completely monotone, so that S(n, t) 
is not a probability distribution. Thus, there is no ran- 
dom process which defines the operational time in such a 
way that the Levy-flight of index ot\ will be transformed 
into a Levy- flight with index ai > a%. The absence of 
an inverse of a TST belonging to a class of subordina- 
tion transformations has a deep physical interpretation: 
a TST is a coarse-graining procedure (see Fig.l): the 
information about the internal steps of the process gets 
lost. One can not anticipate that the transformation in- 
verse to a coarse-graining belongs to the same class as 
the direct one. 
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Fig. 2. The operational time stemming from subordination of 
the two processes depicted in Fig.l. Note that the bar-code- 
like set shows both the intervals of high condensation and the 
long lacunae. 

Note also that the fact that the TET and TST are 
not inverse of each other is mirrored by the fact that 



G 



within the formalism based on the fractional Fokker- 
Planck equations (FFPE), the first one corresponds to 
an additional fractional time derivative in the l.h.s. of 
the FFPE, while the second one i s represented by a frac- 
tional spatial derivative, see Refs. 
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_lj,y,|lj. Note also 
that the noncommutativity mentioned above shows that 
the order of application of these derivatives is fixed and 
cannot be arbitrarily changed. 



VI. THE "PARADOXICAL" DIFFUSION 

A process subordinated to a Levy-CTRW under TET 
(a time transform leading to subdiffusive CTRW) was 
considered in detail in Ref . [ll] . We now know that this 
process subordinates normal random walks under a com- 
bination of TST and TET of different indices (3 and 7. 
The overall behavior of the process is superdiffusive for 
7 < [3 and subdiffusive for 7 > /3. This is easy to un- 
derstand since the scaling considerations show that the 
operational time grows superlinearly with physical time 
in the first case and that the behavior is sublinear in the 
second case. Note that the index fi of the corresponding 
Levy-flight is exactly 2/3, so that this behavior is exactly 
the one obtained in Ref. [[HI. In the case /3 = 7 the op- 
erational time grows linearly with the physical one: Ref. 
jllj suggests that it falls into the diffusion universality 
class. On the other hand this diffusion is a very special 
one: We will call the process subordinated to RW under 
such operational time a paradoxical diffusion. The ran- 
dom process defining an operational time stemming from 
a combination of TST and TET of the same index 7 has 
interesting properties: n typically grows proportional to 
t; on the other hand, neither a well-defined density, nor 
a well-defined mean waiting-time exists. 
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Fig. 3. The PDF of the random walker's positions for para- 
doxical diffusion. The PDFs are plotted as a function of di- 
mcnsionless variable f = x/Q, where Q is the position of the 
upper quartile of the corresponding distribution. The thick 
full line corresponds to RW under the subordination of a TST 
and TET, Eq.Ca); the dashed line corresponds to the inverse 
situation, Eq.(p3fJ. A thin full line represents a Gaussian dis- 
tribution of the same width. 



Let us first discuss the situation mentioned in the 
beginning of the section: a RW subordinated to Levy- 
distributed operational time, driven by a sublinear one. 
The PDF of the corresponding random walks has power- 
law tails, namely, exactly those of a Levy-distribution 
of index 7. On the other hand, the overall width of 
the corresponding curve grows as A ~ \ft. Moreover, 
the whole distribution scales function of dimen- 

sionless displacement £ = x/A: the overall behavior 
is somewhat similar to one found on a polymer chain 
with bridges. The overall form of the function can be 
found using the well-known expression for p(n,u), the 
Laplace-transform of the probability p(n, t) to make ex- 
actly n steps up to time t. Such a process corresponds 
to a directed motion under the same operational time 
as CTRW. For the ordinary renewal process one has 
p(n,u) = i [1- ip(u)]il) n (u), with i(>(u) ~ 1 - u 7 H). 
For u small (t large) this form corresponds to 



p(n, u) 



_1 exp(- 



(26) 



Considering paradoxical diffusion as a process subordi- 
nated to Levy-flights of index 27 under operational time 
given by p(n,t), we get for P(k,u), the Fourier-Laplace 
transform of P(x,t), 



Py(k, li) = I I 

la 



-W 2an p{n,u)dn 



,7-1 



1 27 



U 7 



(27) 



The scaling nature of the distribution is immediately ev- 
ident, the nature of its power-law tails follows from the 
asymptotic analysis for k small: The tail of -P 7 (£) stems 
from those of L(x, 27, 0) and has a power-law asymptotics 
P-y(0 fx £~ 1-27 (7 < 1). Note that such a distribution 
was obtained in Ref. |1T|] as a solution of a fractional dif- 
fusion equation, describing a random process incorporat- 
ing Levy-jumps taking place under sublinear operational 
time. As an example let us consider the distribution 
P 1 / 2 (x,i), i.e. one for 7 = 1/2. This distribution has 
a simple analytical form, which can be obtained by an 
inverse Laplace- Fourier transformation of Eq. (p7j) . The 
inverse Laplace transform of Eq.(|27|) is one given in 3.21 
of Ref. ||and reads: P 1/2 (k, t) 1 = exp(/c 2 i)erfc(|/c| t 1 ' 2 ). 
The inverse (cosine)-Fourier transform of this function is 
given by No. 10.6 of Ref. |24| and reads: 

P\/i{x> t) = ^7r~ 3/2 exp(x 2 /4i)Ei(-x 2 /4t), (28) 

where Ei(x) is the exponential integral, see Eq. 5.1.2 
of Ref. p8| . The corresponding function is a scaling 
function of £ — x/i 1 / 2 ; its behavior for £ large fol- 
lows from asymptotic expansion of — Ei(— x) = Ei(x) = 
x~ 1 e x [l — l/x + ...], so that asymptotically Pi/2(£,) shows 
the £ -2 -like tail, similar to one of Cauchy-distribution. 
For £ — > 00 the distribution P\/2{k,t) shows week (log- 
arithmic) singularity (following from Eq.(5.1.11) of Ref. 
p8[ ) , a sign of strong lacunarity of the corresponding op- 
erational time. The asymptotic analysis of Eq.(|27|) shows 
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that such integrable singularities appear in the center of 
distribution for < 7 < 1/2: the behavior for £ — > 
is given by P 7 (£) ex £ 27_1 , for 7 = 1/2 P 7 (£) diverges 
logarithmically, as we have already seen in Eq.(|2^). 

The distribution Pi/2(£) is plotted in Fig. 3 together 
with the Gaussian distribution (i.e. the distribution 
Pi(£) of the same class, the one corresponding to a nor- 
mal diffusion) and with the distribution stemming from 
the inverse order of application of TET and TST to a 
simple diffusion, which is discussed in detail in the next 
section. All distributions are normalized in such a way 
that their quartiles coincide. Note that the quartiles of 
Pi/a(0 are situated at ±0.841. 



VII. NON-COMMUTATIVITY OF 
TIME-SUBORDINATION 

Applying the transformations other way around, i.e. 
considering a process subordinated to Scher-Montroll 
CTRW under Levy-time, we get a process which is differ- 
ent from one discussed above. Let us start from a simple 
example. 

Let us note that the TET of index 1/2 (correspond- 
ing to an inverse Laplace-transform of the function 
e -' n Vu 1 ' yfu) is given by 



Qi/s(n,t) 



1 



-n 2 /4t 



nt 



(29) 



i.e. corresponds to a part of a Gaussian distribution 
for n > 0, so that n typically grows as t 1 / 2 . The 



corresponding TST is given by a distribution, Eq.(18), 
Ri/2{T,n) = 2 ^ l T 3/i e~" / 4T - The subordination of 
these two processes is described by a function 



S 1/2 (T,t) 



1 



74* _ 



i) 



2 t (T 



2 V ^FT 3 / 2 
-1 



2 / iT dr 



(30) 



which is a probability distribution with the tail decaying 
as (as a tail of a stable distribution of index 1/2) 

and with the square-root singularity at zero. Note that 
this distribution is just a solution of a fractional Liouville 
equation describing directed motion under such an opera- 
tional time, just like Eq.(|28|) is the solution of a fractional 
diffusion equation. This is a process subordinated to a 
Levy one under sublinear time growth. 

We now show that the Q- and R- distributions lead- 
ing to the paradoxical diffusion are not commutative: An 
operatational time resulting from a R*Q transformation 
has a different distribution from one stemming from a 
Q * R-owe. For example, the distribution S\/2(T,t) given 
by Eq.@ is S 1/2 (T,t) = Q* R = J Q(n,t)R(T,n)dn. 
Let us calculate a conjugated disrtibution, S*, 2 (T,t) — 

R*Q = J R(n, t)Q(T, n)dn, one describing a process sub- 
ordinated to a sublinear growth under the operational 



time growing according to a Levy distribution. The dis- 
tribution Sy 2 (T, t) is given by: 



S* 1/2 (T,t) = 



1 e - T2 /4n, _* e - f2 / 4 "rfrc 



20Fn 3 / 2 



2t 1 
~^t 2 + T 2 ' 



(31) 



i.e. corresponds to a positive part of a Cauchy- 
distribution. Note that even such a robust scaling prop- 
erty of a probability distribution as a nature of its power- 
law tail is different from one of the conjugated counter- 
part. 

The plausible scaling consideration here is as follows. 
The distribution Q(T, n) has all moments, so that for n 
large the value of T is well-defined and is of the order 
of n a , a < 1. On the other hand, the distribution of n 
as a function of t is broad and shows a power-law tail 
P{n,t) cx t~ 1 / a {n/t 1 / a )- 1 - a cx in" 1 "". Changing now 
variable from n to T cx n a we get the asymptotics of the 
PDF of Tin a form: P(T, t) cx iT~ 2 , independently on a. 
We note thus that the probability distribution subordi- 
nating a subliner continuous-time directed motion under 
the Levy-distributed operational time of the same index 
has a power-law tail decaying as T -2 , i.e. is similar to a 
Cauchy-distribution. 

The process subordinated to a Gaussian KW under op- 
erational time defined by S* , 2 (T, t) is also not a normal 
diffusion, but represents a marginal situation of a dis- 
tribution whose second moment diverges logarithmically. 
The corresponding PDF shows power-law tails of a x~ 3 
type. This PDF is given by: 



1 



D -x 2 /2T 



2t 



1 



V2^rT 



t 2 + r 2 



dT. 



(32) 



Changing to a new variable £ = x 2 /2T and then in- 
troducing a scaling variable £ = x/^/t we get the PDF 
P(x,t) as a scaling function of £: 



-3/2 



lei 



c 



1/2 p-c 



C 2 +£4/4 



d{. 



(33) 



For £ large the corresponding integral decays as (2/7r)£ -3 . 
Note that Eq.([33"|) can be expressed in terms of Fresnel 
sine- and cosine integrals, S(x) and C(x), so that P(£) 
can be obtained in a closed form: 



p V2(0 - -fe 



s,n 1 2" 



1 - 2S 



1 - 2C 



Id 



(34) 



see Eq.(2.3.7.10) of Ref. f25j| . The corresponding distri- 
bution is also plotted in Fig. 3 as a dashed line. Note 
that the distribution shows a cusp-singularity at £ = 0. 
The value of P(£) in this point is 1/0F = 0.564... The 
quartiles of this distribution are situated at ±0.621. 
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VIII. RELAXATION PHENOMENA UNDER 
TEMPORAL SUBORDINATION 

The fact that the Levy dynamics can follow from a 
temporal subordination is important if one wants to an- 
alyze the possible thermodynamical implications of the 
Levy- flight transport. Imagine an ensemble of thermody- 
namical systems (say, Brownian particles in a harmonic 
potential) which was put out of equilibrium and then let 
relax. As discussed in Sec. 2, such relaxation will lead 
to a stationary state corresponding to a normal equilib- 
rium Boltzmann distribution. Since this distribution is 
time-independent, it would not change under temporal 
subordination, so that the systems with Levy dynam- 
ics may have very ordinary thermodynamical equilibrium 
states and thus be described by normal Gibbs-Boltzmann 
entropy. The non-Boltzmann nature of the equilibrium 
found in Ref. jl9| was connected with the fact that the 
linear response was considered, as proposed by Ref. [^o| , 
an assumption at variance with the findings of Sec. 4. 
Let us now discuss the relaxation to this equilibrium. 

A system slightly outside of the equilibrium can be 
considered as evolving under the influence of the linear 
restoring force. In the operational time of the system 
(marked by the number n of jumps) this relaxation will 
be described by a Fokker-Planck equation. For an over- 
damped particle in a harmonic potential we get, for ex- 
ample: 



On ox \ ox 



(35) 



Note that the values of 7 and D fulfill the Einstein's 
relation, 7 = D/kT. The Green's function of Eq.(||) 
has a form of a Gaussian distribution and reads: 



G(x,n\x ,n ] 
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2vrL>(l - e -27(n-«o)) 

jk(x - e-t^-^xof 
2D(1 - e - 2 7(«-"o)) 



x exp 



, (36) 



see Sec. 5.4 of Ref. p6[ . This equation gives us e.g. the 
PDF at time nina system, in which the particles were all 
situated at x = xq at t — to ■ It is easy to see that the first 
two central moments Mi = (x) and M2 = ((x — (a;)) 2 ) 
relax exponentially to their equilibrium values, so that 



and 



(x(n)} = Xq exp(— t 1 n) 



a 2 (n) = -^-(1 - exp(-2r" 1 n)), 
K7 



(37) 



being a typical pattern of relaxation of a system with 
only one relaxation time r = {k^j)~ 1 . Since all higher 
moments of a Gaussian distribution are the combinations 
of the lower two, they also relax to their equilibrium val- 
ues in a (multi-)cxponential fashion. Let us start from 



6f) and to note that under 



the Fourier-transform of Eq. 
subordination 



P{k, t)= J exp {-ikx e" 7 " - Dk 2 (l - er 2 ~< n )/2 1 ) 

xt- 1/a L(n/t 1/a ,a,-a)dn. (39) 

Let us moreover expand the exponential term in a Tay- 
lor series in k: the coefficients of this series give the mo- 
ments of the corresponding distribution. From Eq.(|39|) 
it follows then that the i-th moment is a combination of 
integrals of the type 



*(*) 



exp 



(^Xn)t-^ a L{n/t^ a ,a,-a)dn (40) 



with A = 7717, < m < i. Using the fact that a Laplace- 
transform of a fully asymmetric Levy-distribution is a 
stretched exponential function, we get: 



= exp(-A(A£ 1/Q ) Q ) = exp(-AA a £). 



(41) 



This means that the exponential relaxation under Levy 
dynamics stays a simple exponential relaxation (only the 
corresponding relaxation time changes). For example, 
the first moment of the distribution (the particle's posi- 
tion) still relaxes exponentially to its equilibrium value 
of zero. On the other hand, the dependence of the relax- 
ation time on the outer parameters (say, temperature) 
entering through the values of 7 and D can change con- 
siderably. Thus, the superdiffusive Levy-flights dynamics 
in the force-free case can coexist with standard thermo- 
dynamics and with very simple relaxation patterns as 
soon as the case of a harmonic force is concerned. 

Let us consider the relaxation in a harmonic potential 
under "paradoxical" diffusion. Here again we can use the 
moment expansion, Eq.(39), and put down the expression 



for the characteristic function of the overall distribution: 



P(k,t) = J exp (-ikx e-~< n - Dk 2 {\ 
x S a (n, t)dn. 



)/27 



(42) 



Note that the moments of the corresponding distribution 
are the combinations of the functions: 



\ exp(-XT)S a {T,t)dT. 
Jo 



(43) 



Note that S a (n, t) is a PDF of a process subordinated to 
a Levy distribution under TET: 



(38) S a (T,t)= dTT-^ a L a (T/r^ a ,a,-a)Q a (T,t)dT (44) 



Thus, a Laplace transform of S according to its outer 
time-variable is a stretched-exponential, so that 



d>(t) 



/ p(T 7 t)exp{-A\ a T)dT 
Jo 



(45) 
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Let us take a Laplace-transform of this expression. Using 
Eq.(^6|) we get: 



$(tt) 



exp(— ru a ) exp(— A\ a r)dr 



(46) 



For small u (long times) this corresponds to a power- 
law decay of 4>(f) of a form <5(t) cx t~ a for t >> A -1 . 
Thus, the relaxation in the case of paradoxical diffusion 
resembles those in normal CTRW and is dominated by 
large lacunae. In the case when the processes are subor- 
dinated other way around, i.e. according to 5*(T, t), the 
decay at longer times follows the universal i _1 -law: for 
example for a = 1/2 we get: 



2/ f°° 1 



2A 



sin(At)ci(Ai) — cos(Ai)si(Ai)] 



(47) 



see Eq.(2.3.7.11) of Ref. |25| (here the integral sine- 
and cosine-functions, si(x) — — J ^^-dx and ci(x) = 
— ^f^dx, are used). For Xt ^> 1 we get: 



(48) 



which asymptotic behavior is universal for all Levy- 
driven CTRWs of the same index. 



IX. CONCLUSIONS 

A broad range of physical processes can be described 
as processes subordinated to a random walk under some 
operational time. In particular, such subordination leads 
to anomalous transport properties, the well-known ex- 
ample being the Scher-Montroll continuous-time random 
walks, a process in which the operational time (given by 
the number of steps) is sublincar in the physical time t. 
Here we have considered the processes subordinated to 
a diffusive process under operational time governed by 
a Levy-distribution with index < a < 1, namely the 
operational time superlinear in physical one. We have 
shown that in the absence of outer forces this subordina- 
tion leads exactly to Levy-flights. The response of such 
a system to weak outer force is strongly nonlinear. In- 
terestingly enough the relaxation patterns in such sys- 
tems are simpler than expected. Thus, we show that the 
behavior in the presence of a weak harmonic force cor- 
responds to a simple exponential relaxation to a normal 
Boltzmann distribution. The combination of super- and 
sublinear operational times (i.e. Levy-flights under sub- 
linear operational time or the Scher-Montroll CTRW un- 
der Levy-time) correspond to the " paradoxical" diffusion, 
a random process which in a force-free case leads to the 
probability-distributions of the particle's displacements, 



which show the power-law tails and lack the second mo- 
ment. The width of the distribution, on the other hand, 
grows proportionally to the square-root of time, showing 
a typically diffusive behavior. Some physical implications 
of these findings have been discussed. 
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